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ABSTRACT 

We have made hydro/N-body simulations with and without star formation to shed some light 
on the conditions under which a central kinematically cold stellar component (characterized 
by a velocity dispersion drop or cr-drop) could be created in a hot medium (e.g. a bulge) and 
survive enough time to be observed. We found that the timescale for a cr-drop formation could 
be short (less than 500 Myr) whereas its lifetime could be long (more than 1 Gyr) provided 
that the central region is continuously or regularly fed by fresh gas which leads to a continuous 
star formation activity. Star formation in the central region, even at a low rate as 1 Mq yr^^, is 
mandatory to sustain a permanent a-drop by replacing heated particles by new low-cr ones. We 
moreover show that as soon as star formation is switched off, the cr-drop begins to disappear. 
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1 INTRODUCTION 

The discovery of a sub-kpc central drop in stellar velocity disper- 
sion (ff-drop in short) in three Seyfert galaxies with double bars 
jEmsellemetai]l200l L hereafter Paper I) has brought this phe- 
nomenon to the foreground. In fact, a a- drop was previo usly ob- 
served in a single object (NGC 6503) by Bottema 1 198 ^. To ex- 
plain t he low velocity dispersion in the centre, Bottem a & GerritsenI 
(l993) had to built up a model with a distinct dynamically cold and 
compact isothermal nucleus in addition to the large scale stellar disc 
with the same mass-to-light ratio. Others examples can be found in 
a few papers ( Heraudeau & Simien 1998; Heraud eau etai]ll999t 
Ijarvis et al. 1988 for instance) but have been most probably over- 
looked because of the low signal-to-noise ratio of these observa- 
tions. More recent observations of larger sample and better signal- 
to-no ise ratio have confirmed I Marauez et al. 2003, Emsellem et al. 
I2OO4I Emsellem et al. in preparation) the presence of a-drop in sev- 
eral single-barred galaxies. However, it is still premature to attempt 
to make statistics on the frequency of a-drops in galaxies since a 
complete and systematic survey is not yet available. More challeng- 
ing and not yet explained is the observation of a-drops in elliptical 
galaxies (e.g. NGC 1404, Gra ham et alJ ll998). The occurrence of 
a-drops in such objects should deserve a detailed study. 

IWozniak et alj (2003) (hereafter Paper II) found that the stel- 
lar aios drop is the kinematical signature of stars bom from a dy- 
namically cold gaseous component. The origin of the drop results 
from the conspiracy of two factors. First, the new stars are con- 
centrated toward the centre, and since they are just bom from the 
low-dispersion gas component, their dispersion is much lower than 
for the old stars component. But second, the effect is amplified by 
the fact that the gas dispersion is also dropping toward the centre 
(and therefore the new stellar component too). This is due to the 
strong accumulation of gas toward the centre, in a nuclear disk. 



where dissipation is stronger than elsewhere, and where the gas is 
therefore cooling down efficiently. 

The existence of flattene d and ra pidly rot ating nuclear stel- 
lar discs has been predicted bvlShlosm an & Be gelmanI il989l) and 
further studied bv lShlosman et air^^9 0). They showed that such 
stellar discs could remain flattened for a long time since two body 
relaxation is a slow process. 

Thus, a few issues remain opened: 

(i) The frequency of a-drop in galaxies depends on the lifetime 
of the phenomenon. In Paper II we were not able to give such an 
estimate. However, even if the frequency of this kinematical effect 
is still under unknown, it can be firmly stated that it is not a pe- 
culiar phenomenon since it is observed in an increasing number of 
'normal' galaxies. We give in Table Q an updated list of a-drop 
observations. It is worth noting that a-drops have been detected 
using a number of absorption lines from the optical to the near in- 
frared (Mg fine triplet at AA 5164, 5173, 5184 A, Call triplet at A 
8542 and 8662 A, ^^CO molecular bands at roughly 2.29^m, etc.) 
which traces various ages of stellar populations. Moreover, a-drops 
cannot be data reduction artifacts since stellar kinematics has been 
retrieved from absorption lines using various techniques. 

(ii) Other mechanisms could be at work to remove some kinetic 
energy from the stellar component. For instance, a massive and con- 
centrated dark halo seems to be able to generate a a-drop in one of 
the simulations of Athanassoula & Misiriotis (2002) (their Fig. 13). 

(iii) Beyond its intrinsic dynamical lifetime, the detectability of 
any a-drop should be also limited by the lifetime of the stellar pop- 
ulations that have been used to trace the kinematics. For instance, 
in Paper I, a-drops has been detected using the ^^CO bandhead at 
2.29 /im. 

In this paper, we will specifically address the problem of a- 
drop lifetime using a new set of N-body experiments. 
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Table 1. List of confirmed or suspect cr-drops in the literature. 



Object 


Morph. Type 


Nucleus 


Reference 


Main abs. lines 


Technique 






Sy 1/2 


M m c 1 1 m f Q 1 / fill S 


Mg b 


pPXF^ 


Nnr 1 nQ7 




Sy 1 


M m c 1 n 1 '/Oil 1 
l-jllldCllCllI CL al. Z.\J\J 1 


i^rn 9 ^iim 


pPXF 


NGC 1138 


SBO 




Simien & Prugniel 2002 


Mp b 


FFM-^ 






ay z 


M m C f*l 1 f*TTI Q 1 / fill 1 


i^rn 9 ^//m 

V- V / Z..J/J;111 


pPXF 


NGC 2639 


SA(r)a 


Sy 1.9 


Marquez et al. 2003 


Call T 


CCF 


NGC 3021 


SAbc 




Heraudeau et al 2003 


Mg b 


FFM 




SBO(s) 




rtgLlClll Cl dl. ZUUj 


Mg b 


Pqq3 


NGC 3fi23 




Liner 


r)p 7ppiiw pt al 7009 

1-/C ZjCCLIW cl al- ^\jyj^ 


Mg b 


FCQ 


NGC 3627 


SBb 


Liner/Sy 2 


Heraudeau & Simien 1998 


Mgb 


FFM 


NGC 4303 


SBbc 


Sy 2 


Heraudeau & Simien 1998 


Mgb 


FFM 


NGC 4579 


SBb 


Liner/Sy 1 .9 


Heraudeau & Simien 1998 


Me b 


FFM 


NGC 4594 


SA(s)a 


Sy 1.9 


Jarvis & Dubath 1988 


Mgb 


CCF* 


NGC 4725 


SBab 


Sy2 


Heraudeau et al. 1999 


Mgb 


FFM 


NGC 4477 


SBO 


Sy2 


Jarvis etal. 1988 


Fe, Mg b 


CCF 


NGC 5728 


SABb(r) 


Sy 2 


Emsellem et al. 2001 


l^CO 2.3Atm 


pPXF 


NGC 6503 


SA(s)cd 


Liner/Hll 


Bottema 1989; Bottema & Gerritsen 1997 


Fe, Mg b 


CCF 


NGC 6814 


SAB(rs)bc 


Sy 1.5 


Marquez et al. 2003 


CallT 


CCF 


NGC 6951 


SAB(rs)bc 


Sy2 


Marquez et al. 2003 


CallT 


CCF 


NGC 7177 


SBb 


Liner/Hll 


Heraudeau et al. 2003 


Mgb 


FFM 


IC 184 


SB(r)a 


Sy 2 


Marquez et al. 2003 


CallT 


CCF 


Circinus 


SAb 


Sy2 


Maiolino et al. 1998 


l^CO 2.3Atm 


CCF 


NGC 3593t 


SA(s)0/a 


Sy2 


Bertola et al. 1996 


5200-6200A 


FCQ 



t edge-on galaxy 

^ Penalized pixel fitting teappellari & EmsellenJl2 004|) 
^ Fourier Fitting Method I Franx. Illingworth & Hec kmanll989t) 
Fourier Con'elation Quotient (Bendei 1990) 

Cross-Correlation Function (e.g. Bottema 1989- and reference therein) 



2 DESCRIPTION OF THE NUMERICAL EXPERIMENTS 

An initial stellar population is set up to reproduce a typical disc 
galaxy. Positions and velocities for 2.5 10^ particles are drawn from 
a superposition of two axisymmetric iMivamoto & Nagail il975h 
discs of mass respectively Mi and M2 (cf. Table|2), of scale lengths 
respectively 1 and 3.5 kpc and common scale height of 0.5 kpc. 
Initial velocity dispersions are computed solving numerically the 
Jeans equations. The initial disc radius is 30 kpc except for run H 
for which the truncation radius is 50 kpc. The gaseous component 
of runs A^*, A400, Agoo and Aiooo is represented by 50 000 particles 
for a total mass of 1.1 10^° Mq distributed in a 6 kpc scalelength 
Miyamoto-Nagai disc. For the collisionless run a"'°°* Mi and M2 
have been proportionally scaled so as to keep the same total mass 
than A°*. 

Evolution is computed with a particle-mesh N-body code 
which includes stars, gas and recipes to simulate star forma- 
tion when necessary. The broad outlines of the code are the fol- 
lowing: the gravitational forces are computed with a particle- 
mesh method using a 3D polar grid with {Nr, N^, Nz) = 
(31,32,64) active cells, the hydrodynamics equations are solved 
using the SPH technique and the star formation process is based 
on Toomre's criterio n for the radial instability of gaseous discs (cf. 
iMichel-Dansac & Wozniak,.2004 for more details). Since we used 
a log-polar grid, we have improved the pre-computation of self- 
forces by subdividing each cell in {nr,n^,nz) = (32,6,6) sub- 
cells. Self-forces are then linearly interpolated before being sub- 
tracted from gravitational forces. 

When star formation is active, the radiative cooling of the 
gas has been computed assuming a solar metallicity. We have also 
checked that our results are independent of the radial resolution in 



the central 100 pc by using Nr, — 40 cells. In Fig.Qwe display 
the region of the bar for various snapshots. More details on the 
global evolution of such kind of simulations, with similar initial 
conditions and star formation on a longer timescale can be found in 
iMichel-Dansac & Wozniafl i2004h for instance. 

Compared to the simulations used in Paper II, those described 
in the present paper are much more accurate since we have used 
5 times more stellar particles and an improved N— body code. We 
will only deal with a case extract from a dozen of such simulations 
of various resolutions and initial setup. However, all our simula- 
tions display a-drops. The generic simulation A^* is thus repre- 
sentative of our d atabase. Another exam ple of simulations with o- 
drops are given in lEmsellem et alj2004l for a study of NGC 1068. 

Apart from the smaller number of stellar particles (1.1 10® 
particles), the setup of the stellar and gas distribution of run H is 
similar to others runs. It includes an additional live dark halo made 
of 2.2 10® live particles distributed in a Plummer sphere of scale- 
length 50 kpc and mass 6.46 10^^ Mq . As our intention is to check 
the robustness of our results against the presence of the massive 
dark halo, a Plummer sphere is a simple but sufficient description 
of the halo. The setup of initial velocities and velocity dispersions 
of disc particles is made consistently with the dark matter distri- 
bution. Since the interaction between the extended live dark mat- 
ter spheroid and the stellar/gaseous disc cannot be computed with 
our particle-mesh code because of the limited vertical range of the 
grid code, run H has been evolve d with GADGET, the tree-based N - 
body+SPH code developed by Sorineel. Yoshida & White (200 ij). 
We used a adaptive time step based on the dynamical time and lim- 
ited by the Courant-Friedrichs-Levy condition. 

H does not exhibit any significant stellar cr-drop up to 5.35 Gyr 
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Figure 1. Evolution of the central ±12 kpc of A' 



Table 2. List of runs 



Model 


Code 


End 

(Myr) 


SF stop 
(Myr) 


Ml 


M2 


Mg 


^nosf 


PM 


2100 


no SF 


0.055 


0.55 




^sf 


PMSPH 


2100 


on 


0.05 


0.5 


0.055 


ssf 
*400 


PMSPH 


643 


400 








71 Sf 

*600 


PMSPH 


801 


600 








7\Sf 

■^1000 


PMSPH 


1612 


1000 








H 


GADGET 


5350 


no SF 









although others kinematical and morphological stellar properties 
are roughly si milar to ha lo-dominated models as those described 
bv lAthanassoula & Misiriotis (2002') (e.g. flat rotation curve, stel- 
lar ring encircling the bar, etc.). The stellar radial velocity disper- 
sion is indeed rather flat, at ft; 180 kms~^ in the central 2 kpc at 
t = 5.35 Gyr. On the contrary, the radial velocity dispersion of the 
gas drops by 20 km s~^ in the central 8 kpc, approximately the di- 
ameter of the stellar bar. Others runs, not displayed here, with more 
massive and/or more extended halos do not develop stellar cr-drops 
either. Since our star formation recipes have not yet been imple- 
mented in GADGET we are not able to check the combined effect 
of a live dark halo and star formation. The effect of a halo on the 
a-drop formation thus deserves a much more in-depth study than 
the present one. This will be discussed in a future paper. 



3 LONG TERM EVOLUTION OF A cr-drop 

To estimate the amplitude of a a-drop (hereafter Act), we have cho- 
sen to express it as the difference between two representative or 
mean values of o: one outside the drop in the kinematically hot 
part of the bulge {o'^'^^) and one inside the drop, i.e. the central 
mean a (cr^'"). Moreover, we will only deal with or for the sake 
of simplicity since the other two components of the stellar veloc- 
ity dispersion (ag, az) also display the same effects (cf. Paper II). 
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Figure 2. Evolution of the cr-drop amplitude estimated as the difference 
between averaged between 450 and 550 pc and the mean inside the 
central 200 pc. The reference run A^* is plotted as a full black line. Star 
formation has been switched off for runs A|qq (blue), A|qq (green) and 
*1000 fr^'') '''' respectively, 400, 600 and 1000 Myr (dashed lines). Vertical 
dotted lines show where star formation has been switched off. Run a""^* is 
also plotted as a full magenta line for reference 

Every value of aji is obviously a mass-weighted average since new 
particles may have different mass than initial ones. 

To determine the value of a^'", we have selected a central 
cylinder of radius 100 pc (i.e. the central 200 pc) at each time. We 
made several trials with height-limited disc instead of a cylinder 
but we found no significant differences apart from a higher noise 
due to lower numbers of particles. The selection typically leads to 
particle numbers between 4000 and 10000 particles in the central 
200 pc. 

After some trials, we decided to compute crJJ'"' in two annuli. 
The first one, between 450 and 550 pc, represents a 'quiet' region, 
outside the drop and the close surrounding region, but still in the 
region that could be considered as the bulge of the galaxy. The sec- 
ond one has been chosen between 250 and 350 pc as being rep- 
resentative of the maximum value of or outside the central drop. 
The difference between o'^'^^ in such a region and a^'" give an 
extreme value of tj-drop. However, the radius of the maximum or 
slightly increases during the evolution. We thus found more practi- 
cal to only use the first annulus (450—550 pc) to analyse into some 
details the typical cr-drop evolution although it slightly underesti- 
mates the real cr^^" and thus the real depth of the a-drop. 

A(T is displayed in Fig.|2| Its evolution could be roughly sep- 
arated into three phases. Boundaries of these phases are approxi- 
mately defined since they should depend on the details of the star 
formation history in the central region. 

3.1 - 300 Myr 

Aa remains approximately constant at 25 kms~^. The comparison 
with the pure N-body run a"°^* shows that such an amplitude is ex- 
pected in case of a low star formation rate or no star formation at 
all. Indeed, it has been shown in Paper II that an artificial a-drop 
of such a small amplitude is expected because of N-body soften- 
ing effects. However, the amplitude of such a numerical a-drop is 
much lower than those observed and should fully diseappear after 
convolution by any realistic point spread function. 

The SFR inside the central 200 pc is displayed in Fig.|3| For 
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whole simulation 
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200 400 600 

Epoch of formation (Myr) 



Figure 3. Evolution of the star formation rate for ran A^* . SFR is displayed 
for the whole simulation (full line), the central 200 pc (dashed hne) and the 
annulus of radii 450-550 pc (dashed-dotted line) 
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Figure 5. Distribution of cru as a function of the epoch of formation for 
particles selected within a radial range [0, 100] pc (top panel) and [450, 
550] pc (bottom panel). Bin size is 25 Myr 



Figure 4. Evolution of the mass inside the central 200 pc for ran A^* . The 
black line is the total mass, the green long-dashed line is the mass of the 
old population, the red short-dashed line is the mass of the new population, 
the histogram-like curve is the mass of the gas. The black dotted line is the 
mass of the new population that is born inside the central 200 pc and is still 
in this region for the time given in abscissa. Vertical dotted lines show the 
three phases discussed in the text 

t < 300 IVIyr the central SFR remains below 1 Mq yr~^while 
most of the new particles are created in the large scale disc. Thus, 
att = 300 Myr, the mass of the new population roughly amounts to 
that of the gaseous component but remains an order of magnitude 
lower than the old population (cf. Fig.|3. 

3.2 300 - 800 IVIyr 

The (T-drop increase is now very significant so that Act reaches 
« 55 km s~^. This is the a-drop growth phase which corresponds 
to the burst of star formation in the central region. However, such 
phase could be due either to the decrease of ctJ^^'" or to an increase 

Ot CT^ . 

To disentangle both effects, we display in Fig. |5| the distri- 
bution of cth as a function of the epoch of particle formation for 
several selected times, up to t — 800 IVIyr. The two regions of 
reference are separately displayed. The old population is included 



in the first bin (0-25 IVIyr). For each curve, the last bin represents 
the mean ct_r of the newest particle bom during 25 IVIyr before the 
corresponding time. 

The an distribution in the central 200 pc clearly shows an in- 
crease with time between t = 200 and t = 600 Myr, for all genera- 
tion of particles. The newest particles have a lower ct_r than the old 
population as expected. At each time, the shape of the histogram is 
roughly a linear relationship between ur and the epoch of particle 
formation. In first approximation, the rate of heating is indepen- 
dent of the epoch of formation, i.e. the age of particles. This heat- 
ing process stops around t — 600 Myr. The ct_r distribution then 
does not evolve any more for particles born between t = 100 and 
t = 600 Myr. A number of particles bom for 50 < t < 100 Myr 
gets into the central region from the outside; these particles have a 
very high ct_r, greater than the old population, but their total mass 
is not high enough to increase ctJ^'". 

Another way to represent ctJ^'" evolution is to display sepa- 
rately the evolution of an for the two particle populations of run 
A''* (a'^'^ for the initial or 'old' particles and ct^^.^"* for the 'new' 
particles, i.e. for particles born during the run) and for all the par- 
ticles (Fig. |6j. The control mn a°°^* is also plotted to show the 
heating of collisionless particles for t > 250 Myr. Such heating is 
mainly due the development of the stellar bar perturbation which 
increases the mass in the central region (cf. Paper II). Run a"°°* 
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Figure 6. Evolution of crj; for particles selected within a radial range [0, 
100] pc (top panel) and [450, 550] pc (bottom panel). Symbols represent 
A°*. Opened squares are the velocity dispersion of the initial population 
full squares are the total cr^ and full circles are the velocity disper- 
sion of the new particles (p'^^ ) created by star formation during the simu- 
lation. Runs A|qq (blue dashed lines), A|qq (green dotted-dashed lines) and 
A|qqq (red dotted lines) can only differ from A^* when star formation is 
switched off after respectively t = 400, t = 600 and t = 1000 Myr (times 
are marked by dotted vertical lines), an for the pure N— body run a"'°°* is 
plotted as a full line 



cannot develop a cr-drop. However, this run allows us to estimate 
the increase of gr due to the secular heating of the central region. 

Fig. |5| (top panel) shows that crj^'" clearly increases until 
t ~ 475 Myr and then slightly decreases to reach a plateau around 
(t55?'"~ 140 km s~^ . The stellar mass in the central 200 pc is mainly 
made out of new particles, especially for t > 400 Myr (Fig|4}. (t^'" 
variations are thus strongly correlated to ctj^'^^ ones. The population 
of new particles responsible for this plateau (those that are bom at 
600 < t < 800 Myr) is made out in part of particles born in- 
side the central 200 pc and in part of low gr particles that have 
migrated from the outside. The fraction of new particles which is 
born and stays in the central region vary from 25 to 10% between 
t = 600 and t = 800 Myr. Thus, the migration of low gr particles 
from outside the central 200 pc takes over from particles born in- 
side the central region when the local SFR strongly decreases after 
t = 475 Myr (cf. Fig.|3)- This dynamical effect is obviously due to 
the strong growth of the central total mass which sinks the potential 
well. 
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The evolution of gr^^ is not driven by the same causes. In- 



Figure 7. As Fig.|5]but for t > 800 Myr. The t = 800 Myr curve is the 
same than for Fig.lsl 



deed, the gr distribution in the 450-550 pc annulus (Fig.|5| bottom 
panel) does not show any saturation of the heating process. On the 
contrary to the central region, there is no clear relationship between 
gr and the epoch of formation. This is obviously due to the fact that 
new particles that are born in the annulus region cannot stay there. 
Thus, both Gr'^ and (J^r"" increase. As a consequence, (Tr^^ con- 
tinuously grows (Fig.|6| bottom panel), its values being moreover 
greatly dominated by gr'^ since the old population represents most 
of the mass in the aimulus. 

The (T-drop growth is thus due in great part to the heating of 
the surrounding region, since (tJJ'" reaches a plateau. 



3.3 > 800 Myr 

During this phase, the cr-drop development seems to saturate. In- 
deed, Act gains only 10 kms^^ in roughly 700 Myr. In the cen- 
tral 200 pc, after t = 800 Myr, the slope of the relationship be- 
tween aR and the epoch of formation changes for particles bom 
after t ~ 500 Myr (cf. Fig0. During this phase, the heating is less 
efficient so that particles bom after t ~ 500 Myr keep a low gr. 
Moreover, most of the mass in the central region is due to particles 
born between 250 and 600 Myr, i.e. during the burst of star forma- 
tion (Fig.|3}. Thus, (tJJ'" does not evolve (cf. Fig|6l, as crj^^", even 
if cr^''' continues to slightly increase. 

The 450-550 pc annulus being dominated by the old pop- 
is driven by the behaviour of gr'^ which slightly 
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increase by a few kms~^. It is noteworthy that crjlp™ is greater 
than crjj'^ during this phase. This effect is mainly due to parti- 
cles bom between 50 and 625 Myr which cross this region. Finally, 
since ag"" does not evolve, the small a-drop fluctuations are again 
mainly due to the weak heating of the surrounding region. 

For t > 1300 Myr, Aa fluctuates around a mean value of 
65 km s~^, with no significant sign of increase or decrease. Thus, a 
(T-drop is not a transient feature. 



4 EFFECT OF SWITCHING STAR FORMATION OFF 

To throw some light on the matter of cr-drop lifetime, we have 
switched off the star formation process at three different times 
of run A^* and pursued the runs (runs are named A%q, Aggo and 
A|qoo) for a while. For A|qq, star formation has been stopped at 
t = 400 Myr, when the total star formation rate (SFR) is maxi- 
mum (roughly 20 Mq yr~^, cf. Fig. |3}. The simulation was then 
proceeded until t ~ 650 Myr. For A|oo> star formation has been 
stopped in the middle of the cr-drop growing phase (t — 600 Myr) 
whereas it has been switched off att — 1000 Myr for Aiqqo, when 
the (j-drop is well developed and closed to be stable. 

The first obvious effect of switching off the star formation is 
the decrease of Act for all the three runs (cf. Fig.|2}. For run A|oo, 
Aa becomes even negative for t > 600 Myr which means that 
the central region becomes hotter than the rest of the galaxy. The 
heating timescale is very short. 

As in Sect.|3| we need to disentangle an variations for both 
representative regions and thus come back to Fig. |6| where (tJJ'" 
and a^fi^^ of A|oo, Agoo and Afooo compared to those of A^*, as 
well as the contribution of cr^''* and a'^" to erg'" and erg'"'. 

For t > 400 Myr, the mass in the central 200 pc is domi- 
nated by the new population when star formation is active, so that 
crg'° closely follows a'^^ . Indeed, the mass of the new popula- 
tion continuously increases when star formation is active (Fig.|4} 
so as to equals the mass of the old population at t ~ 440 Myr. 
Then, the new population dominates the central region. Obviously, 
this is no more the case when star formation is switched off. In- 
deed, between t = 400 and t = 600 Myr, erg'" of Afoo in- 
creases by « 55 km s~^ . The rate at which crg^" increases is higher 
(Ri 0.8 km s"'^ Myr"'^) during the first ^ 50 Myr than the follow- 
ing 150 Myr (~ 0.1 kms^'^ Myr^'^, very close to the rate of run 
^nosf -J ^ marginal effect is the small heating of the old population 
(~ 10 kms"'^ at t = 600 Myr). The heating of the new popula- 
tion is also very efficient in the 450-550 pc annulus but, since this 
region is dominated by the old population, the effect on a'^'^^ is 
limited to ~ 25 kms^'^. Thus the heating of the central region is 
mainly responsible to the cr-drop disappearance. This heating has 
mainly two causes: 

(i) the rapid gas accumulation in the central region strongly in- 
creases the central mass and makes local perturbations of the grav- 
itational potential; 

(ii) low-CT particles that escape from the central region are not 
replaced by new ones. Indeed, the new population is made only 
in part of particles which are bom and remain inside the 200 pc. 
We show in Fig. |4|how evolves the mass of this sub-population. 
It clearly decreases which means that a fraction of low-cr parti- 
cles could be heated outside the central region, during their escape. 
When star formation is active, they are replaced at each timestep by 
the last generation of particles. This is no more the case for A|qq, 
^600 and Af Qoo . 



For A|oo and AiqoO' 'he decrease of Acr has a longer timescale 
than for aIqo even though the same mechanism than for aIqo comes 
into play when the star formation is switched off. The amplitude of 
the decrease of Acr is therefore smaller for Agpg than for run A%q, 
and even smaller for Aiooo ■ Since the central crg^" is roughly con- 
stant for 600 < t < 1000 Myr for A''* (cf. Fig.|S|, one cannot 
invoke the effect of a colder population which needs more energy 
to be heated. In fact, the heating of the central region is much more 
difficult att = 600 or t = 1000 Myr essentially because the mass 
of the new low-cr population is greater than the old one (Fig.|4j- 
Another cause is the bar itself. For t > 475 Myr, the gravitational 
perturbations made by the formation of the stellar bar have van- 
ished. Finally, there is also less gas accumulation in the centre than 
for aIoo since star formation has consumed a significant fraction of 
the available reservoir in the central region. 



5 CONCLUSIONS 

We confirm the result found in Paper II that the appearance of a 
stellar cr-drop in hydro+N— body simulations is unavoidable on the 
condition that star formation is take into account. Indeed, low-cr gas 
accumulates into a central disc because of the gravitational torques 
induced by the stellar bar. It gives raise to a new stellar population 
with the same velocity dispersion. Therefore, the cr-drop appears 
quickly after the beginning of the bar formation, i.e. f» 300 Myr 
for our simulation A^* . It takes roughly 1 Gyr to reach the maxi- 
mum amplitude. Then, the amplitude of the cr-drop remains con- 
stant ( 65 kms^^) until the end of the simulation (2.1 Gyr). 
Timescales could however strongly depend on the star formation 
rate in the central region. 

A cr-drop is thus not a transient feature. However, all stellar 
populations, even low-cr particles, are heated, probably by fluc- 
tuations in the central region of th e gravitational potential fiel d 
as those described by, for instance, jCarlberg & SellwoodlT98^ . 
or stoc hastic heating made by res onances between several density 
waves (Minchev & Guillen 2005). However, in this paper, we did 
not investigate into details what could be the heating mechanisms 
since this deserves a dedicated study. The cr-drop persistence is thus 
not due to some improbable cooling mechanism. Our experiments 
of switching off the star formation recipes at various epoch prove 
the relationship between this long lifetime and the refurbishment of 
the new stellar population of low-cr. A star formation rate as low as 
1 Mq yr"'^ in the central 200 pc is enough to sustain a permanent 
cr-drop by replacing heated particles by new low-cr ones. Such a 
rate is not uncommon in barred galaxies provided that the central 
region is regularly fed by fresh gas. Moreover, when the star forma- 
tion rate is too low, it will take several hundred Myr for the cr-drop 
to disappear. This could explain why cr-drops are still observable 
in early-type galaxies with no sign of nuclear star formation (e.g. 
NGC4477). 

The visibility of a cr-drop as well as its opposition to heating 
mechanisms, strongly depends on the mass ratio between the old 
'hot' population and the 'new' cold one. For A^' the central star 
formation rate reaches 5 Mq yr~^ during the bar formation phase 
so as to create a new population which becomes more massive than 
the initial one alt ^ 440 Myr, thus on a rather short timescale. cr- 
drops should thus be rather common in barred galaxies in the local 
universe. A systematic survey of the central kinematics in barred 
galaxies could confirm these results. 
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